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Abstract 


Weber, J. F., ££ Apestegui, L. (marzo-abril, 2016). Relaciones 
entre parámetros de los modelos de infiltración de Kostiakov 
y Lewis-Kostiakov - Córdoba, Argentina. Tecnología y Ciencias 
del Agua, 7(2), 115-132. 


Se presentan ciertas relaciones halladas entre los parámetros 
del modelo de infiltración de Lewis-Kostiakov, calibrados 
a partir de los resultados obtenidos en una campaña de 
medición in situ en la ciudad de Córdoba, Argentina. Como 
instrumento se utilizó un microsimulador de lluvia portátil 
diseñado y construido por este grupo de trabajo, como 
alternativa superadora al infiltómetro de doble anillo. 
Este equipo, que es completamente desarmable y define 
una parcela de ensayo de 1 m?, permite generar lluvias de 
intensidades comprendidas entre los 65 y 120 mm/h. La 
selección de los sitios de ensayo respondió al tipo y uso 
del suelo. A su vez, se consideraron distintas condiciones 
de humedad inicial. La campaña de medición se extendió 
durante siete meses, realizando entre dos y tres ensayos 
por semana. Los valores medidos fueron digitalizados 
y procesados, construyéndose las curvas de infiltración 
acumulada y tasa de infiltración en función del tiempo. 
A partir de estos valores, se ajustaron los parámetros del 
modelo de infiltración de Kostiakov y de Lewis-Kostiakov 
(Mecenzev), que en sí mismos representan un aporte original 
al conocimiento de la hidrología urbana en esta ciudad. Al 
graficar diversas relaciones entre parámetros, se evidencia la 
existencia de una función empírica, que permitiría disminuir 
a dos el número de parámetros del modelo de Mecenzev, 
con una aceptable pérdida de calidad en las predicciones. Se 
considera que esta relación permitirá mejorar la elección de 
parámetros de infiltración en tareas de diseño hidrológico 
urbano en Córdoba, Argentina. 


Palabras clave: simulador de lluvia, mediciones hidrológicas, 
hidrología urbana, modelo de Lewis-Kostiakov, infiltración. 


Resumen 


Weber, J. E, € Apestegui, L. (March-April, 2016). Relation between 
the Parameters of the Kostiakov and Lewis-Kostiakov Infiltration 
Models, Cordoba, Argentina. Water Technology and Sciences (in 
Spanish), 7(2), 115-132. 


Certain relationships between the parameters of the Lewis- 
Kostiakov infiltration model were found, calibrated based on the 
results obtained from a field measurement campaign in the city of 
Cordoba, Argentina. A portable rainfall micro-simulator was used 
as an instrument, designed and built by this workgroup as a better 
alternative to the double-ring infiltrometer. This equipment, which 
can be completely disassembled, defines a 1 m2 test plot, making it 
possible to generate rainfall intensities between 65 and 120 mm]/h. 
The test sites were selected based on soil type and use. Different initial 
moisture conditions were considered. The measuring campaign 
was conducted over six months with between two and three tests 
per week. The values measured were digitalized and processed and 
accumulated infiltration curves and infiltration rates over time were 
generated. Based on these values, the parameters of the Kostiakov 
and the Lewis-Kostiakov (Mecenzev) infiltration models were fitted, 
which in themselves represent an original contribution to knowledge 
about the urban hydrology in this city. The graphing of the different 
relationships between the parameters indicated the existence of an 
empirical function which would enable reducing the number of 
parameters of the Mecenzev model to two, with an acceptable loss 
in the quality of the predictions. This relationship will improve 
the selection of infiltration parameters for the design of the urban 
hydrology in Cordoba, Argentina. 


Keywords: Rainfall simulator, hydrologic measures, urban 
hydrology, Lewis-Kostiakov model, infiltration. 


Recibido: 04/02/2014 
Aceptado: 22/09/2015 


ISSN 0187-8336 + Tecnología y Ciencias del Agua, vol. VI, núm. 2, marzo-abril de 2016, pp. 115-132 


Weber y Apestegui, Relaciones entre parámetros de los modelos de infiltración de Kostiakov y Lewis-Kostiakov, Córdoba, Argentina 


Introducción 


En trabajos ingenieriles vinculados con el 
diseño y cálculo de obras de desagie pluvial 
urbano, el proyectista se ve en la obligación, a 
la hora de cuantificar los caudales de aporte, 
a considerar dos tipos de superficies en las 
cuencas: impermeables y permeables. Las pri- 
meras sin duda provocan un alto impacto en 
el volumen de crecidas y una disminución del 
tiempo de respuesta de la cuenca (Tucci, 2001). 
En cuanto a las segundas, es posible utilizar 
diversas formulaciones a la hora de cuantificar 
la fracción de la precipitación de escurrimiento 
directo y, por tanto, la fracción infiltrada, entre 
ellas el denominado Método Racional (Chow, 
Maidment, € Mays, 1994), el conocido método 
del CN-SCS (SCS, 1972), muy difundido en su 
aplicación a través de los populares programas 
de computadora HEC-1 (USACE, 1981), y el más 
moderno HEC-HMS (USACE, 2001), el método 
de Horton (1939), difundido en el ámbito de la 
hidrología computacional a través del conocido 
modelo SWMM (entre otros) (USEPA, 1977). Es- 
tos, como muchos otros modelos conceptuales 
del proceso de infiltración (Ravi € Williams, 
1998), intentan simplificar el abordaje riguroso 
del problema, representado por la ecuación de 
Richards (Chow et al., 1994), y por tanto son 
los utilizados más ampliamente en la práctica 
hidrológica. 

Una de las principales dificultades que en- 
cuentra el modelador a la hora de aplicar estos 
modelos de infiltración es la de estimar los 
parámetros de los mismos, tratando de lograr 
una adecuada representación de la realidad. 
Lamentablemente, existen grandes dificultades 
a la hora de disponer de datos medidos de infil- 
tración en áreas permeables urbanas en la Repú- 
blica Argentina. En la ciudad de Córdoba, por 
ejemplo, recién se encara una tarea sistemática 
de producción de información experimental al 
respecto a partir del trabajo de Weber, Urbano, 
Stuyck, Azelart y Martínez (2005), desde el La- 
boratorio de Hidráulica (LH-UTN), dependiente 
del Departamento de Ingeniería Civil de la 
Facultad Regional Córdoba de la Universidad 
Tecnológica Nacional, a través de una campaña 


de obtención de datos experimentales de infil- 
tración en suelos de la ciudad de Córdoba, a 
través de la técnica del infiltrómetro de doble 
anillo. Este procedimiento de ensayo, también 
conocido como método de Miúntz (Custodio 
£z Llamas, 1976), permite obtener mediciones 
directas de infiltración en áreas reducidas (pun- 
tuales). Son conocidas las severas limitaciones 
que presentan las mediciones así determinadas, 
en particular en relación con la falta de represen- 
tatividad del proceso lluvioso en la técnica de 
ensayo, por un lado, y a la fuerte alteración de 
la superficie del terreno debido al hincado de los 
anillos, por el otro. Estas limitaciones conllevan 
a errores experimentales difíciles de cuantificar. 

Con el objetivo de obtener una caracteriza- 
ción experimental más verosímil del proceso 
hidrológico de la infiltración, se propuso utilizar 
como instrumento de ensayo un simulador de 
lluvia portátil desarrollado en el LH-UTN,; se 
visualizan como beneficios inmediatos la mejor 
simulación del proceso lluvia-infiltración y la 
menor alteración de la parcela de ensayo; es 
decir, se intentan subsanar las limitaciones del 
procedimiento del infiltrómetro de doble anillo 
anteriormente utilizado. 


Objetivos 


El objetivo general de esta línea de trabajo es 
la caracterización del proceso de infiltración en 
las áreas permeables de la ciudad de Córdoba, 
Argentina. 

En particular, en este trabajo se propone co- 
mo objetivo la determinación de los parámetros 
de los modelos de Kostiakov y de Lewis-Kos- 
tiakov (Mecenzev) (Ravi €: Williams, 1998) para 
estas áreas, atendiendo a características como el 
uso y tipo de suelo, y con base en las observacio- 
nes experimentales obtenidas mediante el uso 
de un simulador de lluvia; así como explorar 
posibles relaciones entre estos valores, con el 
fin de reducir el número de parámetros y por 
tanto facilitar su selección en tareas de diseño 
hidrológico urbano. Cabe destacar que para 
estos mismos datos experimentales, el grupo de 
trabajo ha determinado los parámetros de otros 
modelos de infiltración, como Horton (Weber, 
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2014a), Green-Ampt (Weber € Apestegui, 2013), 
Philip (Weber, 2015a) y CN-SCS (Weber, 2014b), 
los cuales no son motivo del presente trabajo. 


Metodología 
El aparato simulador de lluvia 


Se describen a continuación las características 
generales del equipo usado en el desarrollo de 
la presente campaña de medición (Weber, Paoli, 
£ Apestegui, 2009, 2010). 

El microsimulador de lluvia tiene la ventaja 
de reproducir con más fidelidad el fenómeno 
de la precipitación. Numerosos autores (Marelli, 
1989) muestran la importancia del impacto de 
las gotas en la microcapa superior del suelo en 
el balance hídrico, además del efecto evidente 
sobre la erosión superficial. Por otra parte, en 
este tipo de dispositivos es posible simular 
una precipitación de intensidad variable en el 
tiempo, incluso una precipitación discontinua. 
La necesidad de reproducir bajo condiciones 
controladas los efectos de la precipitación, y en 
especial el impacto de las gotas de lluvia sobre 
las propiedades físicas del suelo superficial, 
y las consecuencias que ello provoca en la 
infiltración, escorrentía y erosión, condujeron 
desde hace años al desarrollo de diferentes me- 
canismos e instrumentos que tratan de simular 
la lluvia natural (Pla-Sentis, 1981; Rostagno éz 
Garayzar, 1995). El mayor de estos problemas 
recae sobre la reproducción de gotas que se 
asemejen a las de una lluvia natural. 

El formador de gotas utilizado se compone 
de una aguja hipodérmica de 2 pulgadas de lon- 
gitud y 0.8 mm de diámetros interno, introduci- 
da en tubo rígido, el cual se obtiene de cortar su 
propio estuche, al ras de dicha aguja. Se decidió 
colocar un total de 289 formadores, distribuidos 
en 17 filas por 17 columnas, y separados a una 
distancia de 5.30 cm entre ejes, cubriendo de 
esta manera cada formador un área aproximada 
de 28.10 cm? sobre la base de una batea de 96 x 
96 cm de lado. 

El equipo (figura 1) consta básicamente de 
una estructura metálica de caño de acero de sec- 


ción cuadrada, formando un prisma cuadrado 
de 2 m de altura y 1 m de lado. Este prisma se 
forma a través de cuatro columnas y cerrado por 
medio de cortinas rompevientos. Las columnas 
sostienen el sistema de alimentación de agua 
que está compuesto por dos tanques: uno de ali- 
mentación, solamente, y el otro de alimentación 
y regulación de intensidad. La regulación de la 
intensidad se produce a través de un flotante 
(Achutegui, Abreu, £ Páez, 1996). Debajo del 
sistema de alimentación de agua se encuentra 
la batea con los 289 formadores de gotas, que 
producen la lluvia simulada. La base de este 
prisma está constituida por un marco cuadrado 
de hierro de 1 m de lado y 0.10 m de altura, con 
vertederos en sus cuatro lados y canales para 
el transporte del agua escurrida. Este marco se 
hinca en el suelo, constituyendo así las paredes 
límites de la parcela de medición. El peso total 
del equipo armado es de 110.6 kg; el peso del 
fluido puede llegar hasta los 100 kg. 


Figura 1. Simulador de lluvia en campaña. 
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Este equipo dispone de una curva de calibra- 
ción empírica (la cual fue determinada a través 
de un proceso experimental en este laboratorio) 
que relaciona la carga hidráulica disponible 
sobre el extremo de las agujas formadoras de 
gota con la intensidad de lluvia generada. La 
ecuación de regresión obtenida entre estas dos 
variables, con un coeficiente de determinación 
R?= 0.997 fue: 


¡=15.15 08 (1) 


En esta ecuación, ¡es la intensidad de lluvia 
producida (en mm/h) y h es la carga (en cm) 
sobre los formadores. En la figura 2 se presenta 
gráficamente la relación obtenida. 

Del análisis estadístico de los datos de ca- 
libración se pudo concluir que los valores de 
intensidad obtenidos a partir de la ecuación (1) 
tienen una incertidumbre asociada de 3 mm/h, 
lo cual representa un error de entre 3 y 5% de la 
intensidad producida. 


Selección de los sitios 


Al igual que en tareas previas (Weber et al., 
2005), para la definición de los sitios de ensayo 


se consideró que las dos variables a priori más 
influyentes sobre los resultados a obtener eran 
el tipo de suelo y uso del suelo. Con respecto al 
tipo, más que a una clasificación geotécnica, se 
asumió que una clasificación de características 
regionales resultaría con mayor representativi- 
dad hidrológica. Para ello, se fijó como criterio 
la clasificación de ambientes geomorfológicos 
desarrollada por Quintana-Salvat y Barbeito 
(1994). Este trabajo detectó básicamente dos 
ambientes diferenciados en el ejido de la ciudad 
de Córdoba: la planicie loessoide (correspon- 
diente a un suelo del tipo limo arenoso, de gran 
cohesión en condiciones de baja humedad) y la 
terraza de inundación del antiguo cauce del río 
Suquía (primero), la cual fue a su vez clasificada 
por los autores citados en tres subambientes: 
terraza baja, media y alta, según su cercanía al 
actual cauce del río que atraviesa la ciudad. 
En cuanto al uso del suelo, se distinguió 
básicamente en tres tipos fundamentales: patios 
y jardines de viviendas residenciales (espacios 
verdes con bajo tránsito); plazas, parques y 
paseos públicos (espacios verdes con medio 
tránsito); y calles no pavimentadas (calles de 
tierra). El principal criterio de selección se basó 
en la intensidad del tránsito tanto peatonal co- 
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Figura 2. Curva de calibración del simulador de lluvia. 
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mo vehicular debido a su impacto directo en la 
capacidad de infiltración del suelo. Además, se 
agregó como punto de interés especial el lecho 
de una laguna de detención en la zona norte de 
la ciudad. 

En la figura 3 se muestra el plano geomor- 
fológico de la ciudad de Córdoba, con la indi- 
cación de los sitios seleccionados (clasificados 
según el uso de suelo) y los diversos ambientes 
geomorfológicos descritos (Weber, Apestegui, $2 
Baldazar, 2011). Puede observarse la amplitud 
geográfica de la región involucrada. El cuadro 
1 muestra los sitios de ensayo seleccionados 
clasificados según el uso del suelo y el ambiente 
geomorfológico. Como puede verse, no fue po- 
sible conseguir terreno disponible para realizar 
el ensayo correspondiente a espacios verdes de 
bajo tránsito en la terraza media del río Suquía 
debido a que este ambiente cubre un área rela- 
tivamente pequeña y fuertemente urbanizada 
de la ciudad. En total se ensayaron 12 sitios a lo 
largo del ejido urbano de la ciudad. 


Aplicación del simulador de lluvia 


En cada uno de los sitios seleccionados se rea- 
lizaron tres ensayos con el simulador de lluvia 
bajo intensidades relativamente constantes (del 
orden de los 70 mm/h) con diferentes condi- 
ciones iniciales de humedad, las cuales fueron 
obtenidas bien por las condiciones climáticas 
imperantes o porque los terrenos fueron hu- 
medecidos (regados) previamente a partir de 
condiciones secas. Además, en la mayoría de los 
sitios se realizó de modo simultáneo el ensayo 
clásico de infiltración mediante el infiltrómetro 
de doble anillo, como control y contraste de las 
mediciones efectuadas. La campaña de medi- 
ción se extendió durante siete meses e implicó 
el trabajo de tres personas en forma simultánea, 
realizando entre dos y tres ensayos por semana. 

La duración de los ensayos fue variable, 
comprendida en general entre 1.5 y 2 h, depen- 
diendo del comportamiento observado in situ 
del proceso de infiltración, en cuanto a la veloci- 
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Figura 3. Ubicación de sitios de ensayo. 
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Cuadro 1. Sitios de ensayo. 


Tipo de suelo 


Uso del suelo 


(ambiente Espacios verdes Espacios verdes Calles no pavimentadas 
geomorfológico) (bajo tránsito) (medio tránsito) (alto tránsito) 

Terraza alta IT'S Villada ITS Villada ITS Villada 

Terraza media | ==-- Parque General Paz Justo Liebig 5940-B” Villa Belgrano 
Terraza baja Torres de B” Junior Parque de las Naciones Carlos Gauss 4619-B” Villa Belgrano 

. Suá j -B* M. Pp TN 
Planicie loessoide A o ae Hugo Miatello 4600-B” Poeta Lugones 
de Sobremonte Laguna de detención 


dad con que se alcanzó la condición de régimen. 
El número y la frecuencia de las observaciones 
también fue dependiente de las condiciones 
locales encontradas, observándose una menor 
frecuencia de medición en el simulador (y, por 
tanto, una mayor frecuencia en el infiltrómetro) 
en los suelos de mayor capacidad de infiltración. 
Para facilitar la tarea de medición, se seleccionó 
en cada caso el volumen de medición, el cual 
varió desde 50 hasta 500 crm?. 


Determinación del contenido de humedad en 
laboratorio 


Antes del inicio de cada ensayo se extrajeron 
muestras de suelo para determinar el conteni- 
do de humedad antecedente. El contenido de 
humedad fue determinado gravimétricamente, 
a través de la extracción de muestras de suelo 
a 10 cm de profundidad, secándolas en estufa 
a 110 *C y pesando la muestra seca. La hume- 
dad gravimétrica w se determina a través de la 
expresión (2): 


=M (2) 


donde m, es la masa de agua en la muestra de 
suelo, que resulta de la diferencia entre la masa 
de la muestra húmeda menos la masa de la 
muestra seca; m, es la masa de la muestra seca. 
Para relacionar la humedad gravimétrica w (de 
uso principal en geotecnia) con la humedad 
volumétrica O (más frecuentemente usada en 
hidrología subterránea) se utilizó la expresión 
(3) (Tucci, 1993): 
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En esta expresión, d, es la densidad del 
agua (asumida 1 g/cm?) y d, es la densidad 
(peso unitario) del suelo. Este valor se estimó a 
partir de las observaciones in situ realizadas en 
tareas anteriores en cada uno de los ambientes 
geomorfológicos descriptos (Weber et al., 2005). 


Digitalización y proceso de datos 


Los datos obtenidos en la campaña descrita fue- 
ron digitalizados y procesados mediante el uso 
de planillas de cálculo, y después exportados en 
formato de texto plano para su procesamiento a 
través de los códigos desarrollados en específico 
para este trabajo. Cabe destacar que el procesa- 
miento implica también la suma de las series 
temporales de escurrimiento parcial en cada 
uno de los vasos recolectores, para obtener un 
único registro del escurrimiento superficial por 
ensayo. Como una medida del volumen de la 
información experimental generada, se puede 
citar que los 34 ensayos efectuados produjeron 
6 259 pares de valores tiempo-infiltración (un 
promedio de 184 observaciones por ensayo). 


El marco teórico 
El modelo de Kostiakov 


Kostiakov, en 1932 (Ravi € Williams, 1998), 
propuso las siguientes ecuaciones de carácter 
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empírico para estimar la tasa de infiltración y 
la infiltración acumulada: 


F(t)= Ke" (4) 
f (+) =Kat" (5) 


donde K y a son los parámetros del modelo, que 
cumplen las restricciones K> 0 y 0 < a< 1. Es 
de notar que la ecuación (5) tiende a 0 cuando t 
tiende a infinito, por lo que el modelo no suele 
dar buenos resultados para periodos extendidos. 
Además, para las restricciones mencionadas, la 
expresión (5) queda indeterminada para t =0, 
esto es, no es posible definir la tasa inicial de 
infiltración. 


Modificaciones al modelo de Kostiakov 


En razón de lo mencionado, se propusieron 
diversas modificaciones al modelo de Kostiakov, 
en particular el denominado método de Kos- 
tiakov modificado (Al-Azawi, 1985), método de 
Lewis-Kostiakov (Ahuja, Kozak, Andales, $: Ma, 
2007) o método de Mecenzev (Ravi €: Williams, 
1998), que incluye como parámetro adicional la 
tasa base de infiltración f; 


E(t)=f,t+ Ke (6) 
f(t)= f, +Kat” (7) 


En esta formulación, la tasa de infiltración 
tiende a la tasa base L, cuando f tiende a infinito; 
de todos modos sigue siendo indefinido el valor 
de (0). Esto último podría ser significativo en 
la representación de procesos de muy corta 
duración y baja intensidad, que en general no 
tienen aplicación práctica en el diseño hidroló- 
gico urbano. 


Otros modelos 
Otros autores han presentado versiones mejora- 


das del modelo descrito. Clemens (1981, citado 
por Ahuja et al., 2007) describe una mejora del 


modelo de Lewis-Kostiakov mediante la intro- 
ducción de una función escalón sin incrementar 
el número de parámetros del modelo, que ha 
mostrado una mejora sustancial en la predicción 
de la capacidad de infiltración a largo plazo. 


Análisis y definición del modelo a utilizar 


A pesar de las ventajas mencionadas de las 
versiones modificadas, es usual aún encontrar 
aplicaciones de las ecuaciones originales de 
Kostiakov (4), (5) debido principalmente a su 
simplicidad matemática y, a causa de ello, a que 
mediante la aplicación de logaritmos a ambos 
miembros de la expresión (4) se obtenga una 
relación lineal en un papel doble logarítmico: 


InF (+) =InK + alnt (8) 


lo que permite encontrar gráficamente, o a 
través de una simple herramienta de regresión 
lineal, los parámetros del modelo, ya que In K 
es la ordenada al origen y a la pendiente de la 
recta de regresión en papel log-log. 

A diferencia de otros modelos, como el de 
Horton, por ejemplo, donde los parámetros 
tienen una correlación directa con alguna pro- 
piedad de la función [(t), los parámetros K y a 
del modelo de Kostiakov impactan tanto en el 
valor de la tasa inicial de infiltración (entendida 
como el valor de f para un tiempo positivo su- 
ficientemente pequeño) como en la infiltración 
base (entendiendo como tal el valor mínimo que 
alcanza At) en un intervalo de tiempo finito). 
En la figura 4 se muestran diversas curvas de la 
función f(t) para distintos valores de los pará- 
metros K y a. Puede observarse que, a igualdad 
de K, si 4 aumenta, la tasa de infiltración inicial 
disminuye y la tasa base de infiltración sube; y 
que, a igualdad de a, si K se incrementa, tanto 
la tasa inicial como la final de infiltración se 
acrecientan (K representa un factor de escala). 
Para el modelo de Lewis-Kostiakov, el compor- 
tamiento es similar, pero al sumarse la constante 
f, se establece como tasa mínima absoluta ese 
propio valor. 
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Figura 4. Tasa de infiltración según el modelo de Kostiakov para distintos valores de sus parámetros. 


Resultados 


La información experimental así generada fue 
expresada para las variables escurrimiento 
superficial acumulado, infiltración acumulada 
y tasa de infiltración. En las figuras 5 y 6 se 
presentan, a modo de ejemplo, la tasa de escu- 
rrimiento superficial y la infiltración acumulada 
obtenidas para un ensayo en uno de los sitios 
en particular (L. Suárez de Figueroa-B” M. de 
Sobremonte) bajo una intensidad de lluvia de 
66.6 mm/h. 

La variable considerada para el ajuste fue 
la infiltración acumulada, dado que por su 
naturaleza presenta menor dispersión que su 
derivada, la tasa. 


Ajuste de parámetros 


Con base en la información recopilada en cam- 
paña mediante el uso del simulador de lluvia, 
se realizó el ajuste de los parámetros de los 
modelos descritos, buscando la minimización 
de la función objetivo (FO) representada por la 
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suma de las desviaciones cuadráticas entre los 
n valores de infiltración acumulada observados 
(E) y calculados (F*), para cada sitio de ensayo 
y cada condición inicial de humedad: 


F.O.= Y (E -E] (9) 


El proceso se automatizó utilizando un códi- 
go ad hoc desarrollado en GNU Octave (Eaton, 
Bateman, €: Hauberg, 2007), que para el caso del 
modelo de Kostiakov consistió en la aplicación 
de la regresión lineal al conjunto de datos linea- 
lizados por la expresión (8), mientras que para el 
modelo de Lewis-Kostiakov se aplicó el mismo 
procedimiento, pero para la variable reducida 
F' definida como: 


Pont fe (10) 
de forma que el problema del ajuste se reduce al 


de optimación de la variable f,, el que se realizó 
por fuerza bruta con un decimal de precisión. 
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Tasa de escurrimiento (mm/h) 


0.0 0.1 0.2 
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0.6 


Figura 5. Tasa de escurrimiento superficial-Marqués de Sobremonte. 


0.7 


Infiltración acumulada (mm) 


Tiempo (h) 


Figura 6. Infiltración acumulada-Marqués de Sobremonte. 


El método de la fuerza bruta o de búsqueda 
directa (Chapra € Canale, 2006) consiste en 
discretizar el dominio de búsqueda con un 
paso fijo y evaluar la función en la malla así 


generada, detectando de esa forma el mínimo 


buscado. Es claro que la precisión y exactitud 
del método son una función directa del paso o 


intervalo seleccionado. 
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De este modo se obtuvieron 34 juegos de realizado. En el cuadro 2 se presentan los valo- 


parámetros para cada modelo (ya que hubo res ajustados de estos parámetros para cada sitio 
dos mediciones que se consideraron no válidas), de ensayo y para cada medición para el modelo 
junto con el propio valor de la función objetivo, de Kostiakov, y en el cuadro 3 el equivalente 


el cual es una medida de la bondad del ajuste para el modelo de Mecenzev o Lewis-Kostiakov. 


Cuadro 2. Parámetros del modelo de Kostiakov. 


Uso Sitio la medición 2a medición 3a medición 
K=23.01 K= 28.05 K=15.70 
Gauss 4619 a=0.55 a=0.64 a=0.87 
2= 0.995 R?= 0.993 R?= 0.997 
K= 31.68 K=26.50 K= 20.80 
Miatello 4600 a=0.76 a=0.76 a=0.82 
R?2= 0.997 R?=0.997 R?=0.996 
Calles de tierra 
K= 40.06 K = 30.78 K = 25.82 
Liebig 5940 a=0.87 a=0.95 a=0.73 
R?= 0.999 R?= 0.999 R?= 0.999 
K= 20.06 K= 31.05 
TITS Villada a=051l | a=0.75 
R?=0.994 R?=0.996 
K= 50.60 K = 28.97 K= 45.76 
Parque Gral. Paz a=0.87 a=0.72 a=0.85 
R?= 0.998 R?= 0.986 R?=0.996 
K=31.13 K= 30.68 K=31.77 
parque ITS Villada a=0.79 a=0.89 a=0.76 
R?2=0.999 R?= 0.998 R?=0.999 
K= 19.42 K= 29.16 K = 33.44 
Parques Laguna de retención a=0.71 a=0.74 a=0.75 
> 2= 0,997 R?= 0.995 R?= 0.995 
S K= 48.02 K= 33.55 
= Parque de las Naciones a=0.86 | -—- a=0.76 
z R?=1.000 R?= 0.992 
2 K= 33.10 K=32.45 K=37.27 
E UTN a=0.79 a=0.76 a=0.84 
8 R*=1.000 R?=0.999 R?= 0.996 
E K= 43.66 K=27.55 K= 35.67 
E Casa Marqués de Sobremonte a =0.82 a=0.79 a=0.72 
E 
E R?= 1.000 R?=0.998 R?=0.997 
E K= 26.89 K= 34.22 K=40.45 
S Residencias Jardín ITS Villada a=0.81 a=0.71 a=0.76 
z 2= 0.996 R?=0.995 R?= 1.000 
z K=24.20 K=16.65 IESO 
E Torres de Junior a=0.55 a=0.63 a=0.70 
2 R?= 0.998 R?=0.979 R?=0.992 


Te 
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Cuadro 3. Parámetros del modelo de Mecenzev o Lewis-Kostiakov. 


Uso Sitio la medición 2a medición 3* medición 
K= 12.80 K=7.52 K=1.23 
a=0.34 a=0.13 a=0.36 
Gauss 4619 
f,= 10.2 f,=20.4 f,= 14.7 
R?=0.999 R*=1.000 R?= 0.998 
K= 5.25 K= 4.86 K= 1.28 
a=0.19 a=0.26 a=0.05 
Miatello 4600 
f,=26.5 f,=21.7 f,=20.1 
R*= 1.000 R*=1.000 R*= 1.000 
Calles de tierra 
K=5.64 K=2.25 K= 11.20 
haa a=0.42 a=0.81 a=0.50 
Liebig 5940 
f,= 34.4 f,=28.7 f,= 14.6 
R*=1.000 R?= 0.999 R*= 1.000 
K=11.41 K= 5.39 
a=0.29 a=0.19 
IT'S Villada -- 
f,=8.6 f,= 25.7 
R?= 0.999 R?= 1.000 
K= 3.22 K=2.49 K=2.00 
a=0.05 a=-0.20 a=-0.09 
Parque Gral. Paz 
f,= 47.6 f,=26.7 f,= 4.2 
R?= 1.000 2= 0.999 R?= 0.999 
K=8.01 K=0.47 K= 10.95 
a=0.42 a=-0.39 a=0.47 
Parque ITS Villada 
f,=23.1 f,=30.5 f,=20.8 
R?= 1.000 R*=1.000 R?= 0.999 
K= 6.93 K= 5.92 K= 5.16 
a=0.39 a=0.24 a=0.10 
Parques Laguna de Retención 
f,= 12.4 f,=23.4 f,= 28.3 
R?= 0.998 R*=1.000 R?= 1.000 
K= 14.81 K=2.75 
a a=0.63 a=-0.10 
Parque de las Naciones -- 
f,= 33.2 f,=31.1 
R*= 1.000 R*= 1.000 
K= 18.46 K=32.45 K=2.31 
a=0.66 a=0.76 a=-0.18 
UTN 
f,= 14.6 f,=0.0 f,= 34.8 
R?= 1.000 R?=0.999 R?= 0.999 
K=27.91 K=8.24 K=9.66 
a=0.74 a=0.49 a=0.27 
Casa Marqués de Sobremonte 
f,=15.7 f,=19.3 f,=25.9 
R?= 1.000 R?=0.999 R?= 1.000 
K= 26.89 K=9.67 K= 40.45 
a=0.81 a=0.30 a=0.76 
Residencias Jardín IT'S Villada 
f,=0.0 f,=24.5 f,=0.0 
R?= 0.996 R?= 0.998 R?= 1.000 
K= 23.80 K=2.53 K=4.52 
a=0.54 a=-0.07 a=0.11 
Torres de Junior 
f,=0.4 f,=14.3 f,= 19.4 
R?= 0.998 R?= 0.998 R?= 0.999 
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Como puede verse, los valores de R? son 
altos, comprendidos entre 0.9301 y 0.9998 para 
el modelo de Kostiakov, y 0.9575 y 0.9999 para 
Lewis-Kostiakov. De hecho, y como era de espe- 
rar, se obtienen mejores ajustes con este último, 
lo que se ve reflejado en la figura 7, donde se 
presentan los valores de R* obtenidos con el 
modelo de Lewis-Kostiakov, en función de los 
R? obtenidos con Kostiakov; puede verse que el 
conjunto de puntos se ubica por encima de la 
recta identidad. 

Es de notar que en algunas mediciones, el 
parámetro a del modelo de Lewis-Kostiakov 
resultó negativo; esto se explica porque la capa- 
cidad de infiltración creció ligeramente durante 
el ensayo, pudiéndose justificar por una ligera 
disminución en la intensidad de la lluvia gene- 
rada o a particularidades en el sitio de ensayo. 

En las figuras 8 y 9, a título de ejemplo, se 
presenta la infiltración acumulada, resultado de 
una medición en particular (primera medición 
calle Gauss), junto con el ajuste obtenido con 
los modelos de Kostiakov y Lewis-Kostiakov, 
respectivamente. En los cuadros 4 y 5 se presen- 
tan los valores medios (junto con su rango de 


variación) de los parámetros de los modelos de 
Kostiakov y Lewis-Kostiakov, respectivamente, 
en función del uso del suelo. 

En todos los casos, la inclusión del término 
f,t en la expresión de Lewis-Kostiakov de la 
infiltración acumulada (6) se tradujo en una 
reducción de los valores de los parámetros K 
y a, como puede observarse en las figuras 10 
y 11. 

Se analizó una posible relación de depen- 
dencia entre los parámetros de ambos modelos, 
graficando, para el modelo de Kostiakov, K 
versus a, y para el modelo de Lewis-Kostiakov, 
K versus a, K versus f,, f, versus a; observándose 
una posible relación en el caso de K versus a para 
Lewis-Kostiakov (figura 12). En estos gráficos 
también se clasificaron los puntos de acuerdo 
con el uso del suelo, no pudiéndose identificar 
un agrupamiento a partir del uso, aunque se 
aprecia que los valores del parámetro K del 
modelo de Lewis-Kostiakov para calles no pavi- 
mentadas presenta un rango más bajo y acotado 
que para los otros usos del suelo (figura 12). 

En la figura 12 puede apreciarse un valor 
atípico (outlier) en K=2.2488 y a =0.8125. Extra- 


1.000 


0.995 


0.990 


0.958 


R? Mecenzev 


0.980 


0.990 0.995 1.000 


R? Kostiakov 


Figura 7. Coeficiente de determinación R? obtenido con el ajuste del modelo de Mecenzev en función del R? conseguido con el 


ajuste del modelo de Kostiakov. 
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Figura 8. Infiltración acumulada, medida y obtenida por el modelo de Kostiakov-primera medición, calle Gauss. 
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Figura 9. Infiltración acumulada, medida y obtenida por el modelo de Lewis-Kostiakov (Mecenzev)-primera medición, 
calle Gauss. 
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Cuadro 4. Parámetros medios del modelo de Kostiakov según el uso del suelo. Entre paréntesis, el rango de variación. 


Uso K a 
Calles 26.68 (15.70-40.06) 0.75 (0.51-0.95) 
Parques 34.67 (19.42-50.60) 0.79 (0.71-0.89) 
Viviendas 30.35 (16.65-43.66) 0.72 (0.55-0.82) 


Cuadro 5. Parámetros medios del modelo de Lewis-Kostiakov según el uso del suelo. Entre paréntesis, el rango de variación. 
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Uso K a f, (mm/h) 
Calles 6.26 (1.23-12.80) 0.32 (0.05-0.81) 20.51 (8.60-34.40) 
Parques 8.28 (0.47-32.45) 0.20 (-0.39-0.76) 26.48 (0.00-47.60) 
Viviendas 17.07 (2.53-40.45) 0.44 (-0.07-0.81) 13.28 (0.00-25.90) 


K Mecenzev 


K Kostiakov 


Figura 10. Relación entre los parámetros K de los modelos de Kostiakov y de Lewis-Kostiakov (Mecenzev). La línea continua 


corresponde a la función identidad. 


yendo este valor de la serie fue posible ajustar 
una expresión de la forma: 


a=alnK-$ (11) 


con a = 0.261145, f = 0.203251 y R*=0.76586. A 
partir de este resultado, se planteó entonces la 
posibilidad de una modificación al modelo de 
Lewis-Kostiakov, denominado LK-2p (Lewis- 
Kostiakov a dos parámetros) de la forma: 
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Elt)= HRK (12) 


donde f, y K deberían ajustarse en cada medi- 
ción, y a y fB serían constantes en el nivel regio- 
nal. Para explorar esta hipótesis, se desarrolló 
un código ad hoc en Octave que permitió ajustar 
de modo global (es decir, sobre las 34 medicio- 
nes) los valores óptimos de a y f, y localmente 
(esto es, para cada una de las 34 mediciones) los 
valores óptimos de f, y K. La optimación global 
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Figura 11. Relación entre los parámetros a de los modelos de Kostiakov y de Lewis-Kostiakov (Mecenzev). La línea continua 


corresponde a la función identidad. 


O Calles 


O Parques 
X Viviendas 


a Mecenzev 


K Mecenzev 


Figura 12. Relación entre los parámetros a y K del modelo de Lewis-Kostiakov (Mecenzev). 


de a y [3 se realizó por medio de un autómata aproximada central en las ocho direcciones 
celular sencillo, minimizando la suma de los posibles, identificando la posición del mínimo, 
residuos cuadráticos (ecuación (9)). En este reposicionando la malla, y repitiendo el proce- 
procedimiento, un apuntador recorre una malla so hasta que la posición del mínimo no varíe. 
regular bidimensional a partir de una solución La optimación local de los valores de f, y K se 
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realizó aplicando el método de Nelder-Mead 
(Mathews é Fink, 2004). En ambos casos se hizo 
una Optimación no restringida. Como valores 
óptimos se obtuvieron a = 0.2560 y f = 0.2040, 
valores relativamente cercanos a los originales. 


Los valores locales de los parámetros f, y K 
pueden observarse en el cuadro 6. Si bien en 
algunos casos los parámetros así ajustados dan 
lugar a inconsistencias —valores negativos— en 
un valor de K y en tres valores de f,, en la ge- 


Cuadro 6. Parámetros del modelo LK-2p. 


Uso Sitio la medición 2a medición 3a medición 
K=22.78 K= 14.42 K=-0.90 
Gauss 4619 f, =0.40 f,= 18.73 f,=17.29 
R*= 0.997 R?*= 0.998 R*=0.981 
K=5.81 K=4.07 K= 1.00 
Miatello 4600 f, = 25.98 f,= 22.45 f, = 20.00 
tin R?=0.999 R*=1.000 R?*=0.999 
K=3.27 K=2.06 K= 7.73 
Liebig 5940 f, = 36.66 f,= 29.30 f,=17.99 
R?= 0.999 R?= 0.999 R*=1.000 
K=20.29 K= 6.15 
IT'S Villada f,= 0.08. | f, = 25.00 
R?=0.996 R*=1.000 
K=3.75 K=4.55 K=2.40 
Parque Gral. Paz f£ 47.4 f,= 24.80 f,= 43.91 
R*= 1.000 R?*= 0.998 R?=0.999 
K=5.57 K= 28.97 K=7.79 
Parque ITS Villada f,= 25.44 f, = 0.89 f, = 23.90 
R*= 1.000 R*= 1.000 R*=0.999 
K=4.84 K=6.15 K=7.54 
Parques Laguna de Retención f,= 14.42 f,= 23.19 f,=26.01 
pa R?= 0.998 2= 1.000 R*=1.000 
2 K=5.95 K=4.36 
y Parque de las Naciones O f,= 29.59 
£ R?= 1.000 R?=0.999 
2 K=8.53 K=18.99 K= 3.59 
E UTN f, = 24.38 f,=183.12 f, = 33.66 
8 R?= 0.999 2= 0.999 R*=0.999 
E K=12.86 K=5.18 K=15.07 
E Casa Marqués de Sobremonte f, = 30.52 f,=22.18 f,= 20.64 
El 
E R?=1.000 R?= 0.999 R*=1.000 
E K=4.29 K=13.64 K=54.31 
S Residencias Jardín IT'S Villada f, = 22.58 f,=20.61 f,=-13.80 
z R?= 0.993 R?*= 0.998 R*=1.000 
z K= 35.04 K= 3.38 K= 5.57 
E Torres de Junior f,=-10.76 f,=13.51 f,= 18.40 
S R?= 0.998 R*= 0.998 R*=0.999 
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neralidad los resultados son aceptables. Como 
es de esperar, los valores de R? son, en general, 
ligeramente menores para LK-2p frente a Lewis- 
Kostiakov (sólo en tres de los 34 ensayos R?es 
mayor para LK-2p), pero la reducción de tres a 
dos parámetros locales libres podría justificar 
no tener tanta exactitud. 


Conclusiones y recomendaciones 


A partir de la información recopilada in situ 
durante la campaña llevada a cabo en diversos 
sitios de la ciudad de Córdoba ha sido posible 
ajustar los parámetros de los modelos de infiltra- 
ción de Kostiakov y Lewis-Kostiakov (Mecen- 
zev), resultando en la mayoría de los casos este 
último como el de mejor desempeño. Al analizar 
una relación entre los parámetros a y K de este 
último, se observó una tendencia a una relación 
logarítmica, cuyas constantes de carácter global 
fueron ajustadas de manera simultánea con los 
parámetros locales f, y K, dando lugar a un mo- 
delo de dos parámetros denominado LK-2p, que 
presenta algunas inconsistencias, producto de la 
no imposición de restricciones al ajuste; pero, 
en términos generales, se obtienen resultados 
aceptables, con valores de R? superiores a los del 
modelo de Kostiakov y ligeramente inferiores a 
los de Lewis-Kostiakov. 

Si bien los modelos de Kostiakov y Lewis- 
Kostiakov no representan las formulaciones 
algebraicas más sofisticadas para la representa- 
ción del proceso de infiltración, sí han mostrado 
un buen ajuste a los datos experimentales reco- 
lectados, incluso superior al de otros métodos 
de amplio uso, como el del CN-SCS (Weber, 
2015b), y resulta atractiva la idea de aumentar 
la parsimonia a través de una relación empírica 
como la que da origen al denominado modelo 
LK-2p. 

Resta investigar si la relación logarítmica (12) 
es propia de la información local obtenida en las 
mediciones realizadas en la ciudad de Córdoba 
o responde a un carácter más general. 

En cualquier caso, se considera que el con- 
junto de parámetros así ajustados será de interés 


para los proyectistas vinculados con el manejo 
del escurrimiento pluvial urbano en la ciudad. 
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